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Amendments to the Claims 

This listing of claims will replace all prior versions, and listings, of claims in this application. 

The following amendments do not constitute an admission regarding the patentability of the 
amended subject matter and should not be so construed. Amendments to the claims were made for 
purposes of more clearly stating the claimed subject matter and do not add new matter or alter the 
scope of the claims. 

Listing of Claims; 

1 . (Currently amended) A method for constructing a gene network of relationships 
between a plurality of genes of a set of genes and generating a graph representing the gene network 
of relationships between the genes comprising the steps of: 

(a) providing a quantitative time coxjrse data library for tiie a set of genes of an organism, 
said library including expression results based on time course of expression of each gene in said set 
of genes, wherein the expression results comprise gene expression data at a plurality of time points 
for each gene in said set of genes, quantifying an average effect on expression of each gene at each 
time point by each other of said genes at each time point, and quantifying a measure of variability of 
the effect on expression of each gene at each time point on each other of said genes; 

(b) creating a gene expression matrix from said librar y, wherein the gene expression matrix 
comprises matrix data comprising the effects on expression of each gene at each time point by each 
other of said genes at each time point : 

(c) a pplying a Bayesian computational model to the matrix data, wherein said Bavesian 
model comprises minimizing a BNRC/yv„/,«,> criterion: and 

(e) (d) generating the gene network of relationships between said genes and generating a 
graph representing the gene network of relationships between said genes.f «td 

(d) determining if on e or mor e groups of g e n e s is e xpr e ss e d diff e r e ntly from oth e r of said 
groups of g e n e s. 
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2. (Cancelled) 

3. (Currently Amended) The method of claim i\, wherein said step of minimizing a 
BNRCrfy„a;„,c criterion comprises using a non-linear curve fitting method selected from the group 
consisting of polynomial bases, Fourier series, wavelet bases, regression spline bases and B-splines. 

4. (Currently amended) The method of claim 1, wherein said data library is created by 
obtaining the gene expression data at the plurality of time points from a a sing time course study te 
alt e r g e ne e xpr e ssion, studying altering gene expression over time. 

5. (Currently amended) The method of claim 31, wherein said step of minimizing said 
BNRCrfy„amic criterion further comprises selecting a Bayesian mode using a backfitting algorithm. 

6. (Currently amended) The method of claim 21, wherein said step of minimizing a 
BNRCrfj,„aOTic criterion further comprises using Akaike's information criterion. 

7. (Currently amended) The method of claim 31, wherein said step of minimizing a 
B'NR.Cdynamic Criterion further comprising comprises using maximum likelihood estimation. 

8. (Currently amended) The method of claim 1, wherein said plurality of genes are 
associated with a cell cycle of the organism . 

9. (Currently amended) The method of claim 21, wherein said measure of variability is 
variance. 

10. (Currently amended) The method of claim 21, wherein said non- linear curve fitting 
method is a non-parametric method. 

1 1 . (Original) The method of claim 1 0, wherein said non-parametric method for 
minimizing a B'NKCdynamic criterion includes using heterogeneous error variances. 



12. (Currently amended) The method of claim 11, wherein said step of minimizing a 
BNRCrfynaOTic criterion further comprises the steps of: 
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(1) making a score matrix whose (i, j)* element is the "Si^O dynamic score of the graph 

gene,-^ geney; 

(2) implementing one or more of add, r e move and revers e adding the BNRO j iTOTmjr score 
of the effect of parent gene/ on gene^ to a calculation of BNRC^v„y,m,v-, removing the B'NKOHy «nn^ 
score of the effect of parent gene, on gene^ to a calculation of BNRCrfv»/,«,>. and reversing the 
BNRC^//v»„„,> score of the effect of parent gene; on gencj to a calculation of ENRC^ i;;^. which 
provides the smallest BNRCrf^„a;„,c; and 

(3) repeating step 2 until the BNRCflr_y„aOT(c does not reduce further. 

1 3 . (Original) The method of claim 1 1 , wherein said step of minimizing a BNRCrfy„cm«c 
criterion further comprises the step of applying a hill-climbing algorithm to minimize dynamic- 

14. (Currently amended) The method of claim 11, wherein an intensity of the an edge is 
determined using a bootstrap method. 

15. (Ciurently amended) The method of claim 14, wherein said bootstrap method 
comprises the steps of: 

( 1 ) providing a bootstrap gene expression matrix by randomly sampling a number of 
times, with replacement, from the original gene library expression data; 

(2) estimating the gene gen e tic network for genCj and genej; 

(3) repeating steps (1) and (2) 7" times, thereby producing T gene geaetie networks; and 

(4) calculating the bootstrap edge intensity between genci and genej as (ti+t2)/T. 

16. (Withdrawn) A method for elucidating a gene network, comprising the steps of: 

(a) providing a raw data library of time-course gene expression data for a plurality of 
genes of an organism; 

(b) subtractmg background signal mtensities from said raw data library; 
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(c) calculating the relative change in gene expression for each of said plurality of genes; 

(d) analyzing the statistical significance of said relative in gene expression using 
Student's t-test; and 

(e) fitting said changes in gene expression to a linear spline function. 

17. (Withdrawn) The method of claim 16, fiirther comprising the step of removing from 
consideration, those genes whose expression levels are sufficiently low so as to be determined 
predominantly by noise. 

18. (Currently amended) The method of claim 1, wh e r e in said step of grouping 
comprises comprising grouping said genes into one or more equivalence sets. 

1 9. (Withdrawn) A method for estimating a gene network relationship, comprising the 
steps of: 

(1) Making a p x p matrix whose (i j)th element is a BNRC score of the graph genci -» 

genej; 

(2) select a candidate set of parent genci of genej that gives a small BNRC score 

(3) select a computational order of said parent genes; 

(4) repeat the following steps; 

(4. 1) for genes either add a parent gene or delete a parent gene; 

(4.2) recalculate BNRCd3,„„«ic score; 

(4.3) repeat steps 3.1 and 3.2 until a suitable convergence criterion is satisfied; 

(5) permute the computational order of said parent genes in step (3); 

(6) repeat step (4); and 

(7) repeat steps (5) and (6) until B'NR.Cdynamic is minimized. 
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20. (Currently Amended) A method for constructing a gene network model of a system 
containing a network of genes from time course gene expression data, said method compris e s 
comprising: 

applying using a Bayesian computational model to the time course gene expression 
data, wherein a pplying said Bayesian computational model comprises minimizing a 
BNRCrf,,B/,m//- criterio n: and 

generating a graph representing the gene network model resulting from applying the 
Bayesian computational model to the time course gene expression data wherein BNRC, rfi;;»^ 
criterion is minimized . 

2 1 . (Original) The method of claim 20, wherein minimizing the B^Cdynamic criterion 
comprises using a non-linear curve fitting method selected from the group consisting of polynomial 
bases, Fourier series, wavelet bases, regression spline bases and B-splines. 

22. (Currently amended) The method of claim 20, wherein minimizing the BNRCdy„a„,c 
criterion comprises selecting a Bayesian mod e model using a backfitting algorithm. 

23. (Original) The method of claim 20, wherein minimizing the B'i^Cdynamic criterion 
comprises using Akaike's information criterion. 

24. (Original) The method of claim 20, wherein minimizing the B'NS.Cdynamic criterion 
comprises using maximvun likelihood estimation. 

25. (Original) The method of claim 20, wherein minimizing the Bb^KCdymmic criterion 
comprises using a non-linear curve fitting method, wherein the non-linear curve fitting method is a 
non-parametric method. 

26. (Original) The method of claim 25, wherein the non-parametric method includes 
using heterogeneous error variances. 

27. (Currently amended) The method of claim 26, wherein minimizing the BmiCdynamic 
criterion ftirther comprises the steps of: (1) making a score matrix whose (i, j)* element is the 
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BNRC^/y„amic score of the graph gene/ -> gene,; (2) implementing one or more o f add, remov e and 
Fevefse adding the BNRC^rf ,^gg ^ score of the effect of parent gene/ on p:enej to a calculation of 
BNRC^v„„m/r, removing the BNRd score of the effect of parent gene/ on gene^ to a calculation 
of BNRCrfv„.,^,>. and reversing the BNRdrfw„„;r score of the effect of parent gene,- on gene^ to a 
calculation of BNRC/A;»<7»,/r. which provides the smallest BNKCdynamic, and (3) repeating step 2 xmtil 
the Bl^iRCdynamic docs not reduce fiirther. 

28. (Currently amended) The method of claim 26, wherein minimizing the B'i>iR.Cdynamic 
criterion further comprises the step of applying a hill-climbing algorithm to minimize the 

B-m.C^'^ dynamic- 

29. (Currently amended) The method of claim 26, wherein an intensity of fee an edge is 
determined using a bootstrap method. 

30. (Currently amended) The method of claim 29, wherem said bootstrap method 
comprises the steps of: 

(1) providing a bootstrap gene expression matrix by randomly sampling a number of 
times, with replacement, from the original gene library expression data; 



(2) 


estimating the gene genetic network for genei and genej; 




(3) 


repeating steps (1) and (2) Ttimes, thereby producing Jgene get 


etie networks; and 


(4) 


calculating the bootstrap edge intensity between genCj and gene, a 


S(ti+t2)/T. 


31. 


(Cancelled) 




32. 


(Cancelled) 




33. 


(Cancelled) 




34. 


(Cancelled) 
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35. (Withdrawn) A method for identifying a target gene in a system containing a gene 
network, comprising: (a) constructing a first and second gene network model using a Bayesian 
computational model, wherein said Bayesian computational model comprises minimizing a 
BNRCdynamic critcrion, wherein the first gene network model is obtained by analyzing a first gene 
expression profile and the second gene network model is obtained by analyzing a second gene 
expression profile, and wherein the first gene expression profile is obtained fi-om the system at a first 
time point and the second gene expression profile is obtained fi'om the system at a second time point 
after said first time point, and (b) analyzing the first and second gene network model using said 
Bayesian computational model, wherein the time course of gene expression is quantified, and 
wherein a parent gene is identified as the target gene. 

36. (Withdrawn) The method ofclaim 35, wherein the target gene is a parent gene. 

37. (Withdrawn) The method of claim 35, wherein the target gene is a gene downstream 
fi"om a parent gene. 

38. (Withdrawn) A data file containing the identity of one or more target genes obtained 
according to the method of claim 35. 



39. 


(Withdrawn) 


The data file ofclaim 38 in a computer readable form. 


40. 


(Withdrawn) 


The data file of claim 38 accessible fi"om a remote location. 


41. 


(Withdrawn) 


The data file ofclaim 38 accessible fiom an internet web location. 


42. 


(Withdrawn) 


A method of providing a service comprising (1) receiving a data set 



from a party, said data set comprising time course expression data of a group of genes; and (2) 
determining network relationships between genes in said group by minimizing a Bb^Cdymimic 
criterion. 

43. (Withdrawn) The method of claim 42, wherein receiving said data set comprises 
receiving the identity of at least one of said genes. 
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44. (Withdrawn) A method of providing a service comprising receiving an agent from a 
party, and identifying a target gene for the party using the gene network model constructed 
according to the method of claim 35. 
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